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Summary 

Equations for the power absorbed by a glass tube 
have been derived and applied to the particular case 
of a highly transmissive, fused silica tube irradiated 
by concentrated sunlight in space near Earth. Cal- 
culations were made with the assumption that the 
focal volume of the solar concentrator entirely en- 
velopes the tube. Tube parameters studied were out- 
side diameter and wall thickness. Absorption data 
between 155 and 4900 nm were used. Absorption 
beyond 4900 nm was estimated. Volumetric radia- 
tion and the variation of absorptance and emissivity 
with temperature have been used to show that, with 
radiative cooling only, large solar concentrations can 
develop very high temperatures in a fused silica tube. 

Introduction 

NASA has been investigating solar-pumped lasers 
(refs. 1 and 2) and is now studying space-based 
concepts for such lasers (refs. 3 and 4). One proposed 
container for the lasing medium of such a system is 
a fused silica tube that is to be placed at the focus 
of a large reflector capable of concentrating sunlight 
several thousand times at the surface of the tube. 
Since such concentrations of sunlight may produce 
high temperatures in the tube, a means is needed 
to estimate that temperature. This paper provides 
that means in the form of equations that can be 
used to calculate the power absorbed in glass tubes. 
Conversion of the absorbed power to temperature 
is determined for the particular case of a highly 
transmissive form of fused silica. 

Analysis 

In the derivation of equations which follows, ex- 
tensive use has been made of the attenuation fac- 
tor e, which is a function of the absorption coeffi- 
cient and the path length in glass at a particular 
wavelength. This absorption coefficient may be cal- 
culated from spectral data on the index of refraction 
and the spectral transmission of the glass (for nor- 
mally incident light) by use of the following formulas, 
which are derived from multiple reflection, transmis- 
sion, and attenuation phenomena in flat glass sheets 
of thickness d : 


R = 



( 1 ) 


T (1 ~R) 2 e 
1 - e 2 R 2 


( 2 ) 


( 3 ) 


where 


R 

reflectance of interface between 
indices of refraction 1 and n 

a 

absorption coefficient 

T 

intensity transmission 


(The above equations may be found in other litera- 
ture (e.g., ref. 5), but are presented here for conve- 
nience and continuity in the following analysis.) 

Reflection from, transmission through, and atten- 
uation in a glass tube are not as simple as those of a 
flat glass plate. If the image of the Sun at the tube 
is greater than or equal to the outside diameter of 
the tube, then at any point on the tube surface rays 
of sunlight enter the tube at all angles of incidence 
from —90° to 90°, as shown in figure 1. (All rays 
are in a plane perpendicular to the tube axis.) The 
ray paths inside the tube are in three angular re- 
gions: region 1 (fig. 1(a)), where the rays penetrate 
the inner surface; region 2 (fig. 1(b)), where the rays 
are totally internally reflected from the inner surface; 
and region 3 (fig. 1(c)), where the rays never reach 
the inner surface. The reflectance at each surface is, 
in general, different because the angle of incidence is 
different. Absorptance at each entrance angle is dif- 
ferent because the path lengths of the corresponding 
rays differ. An infinite series of reflections is formed 
between the inner and outer surfaces beginning at 
any point where a beam enters the glass. In angular 
region 1, for example, a beam enters the outer sur- 
face once and parts of it enter the inner surface an 
infinite number of times, producing at each entrance 
an infinite, diminishing series of reflections between 
inner and outer surfaces. In angular regions 2 and 3, 
beams do not penetrate the inner surface, but an in- 
finite series of reflections is generated by the beam 
entering the outer surface. 

The key to analyzing this complex model lies in 
analyzing one series of reflections between inner and 
outer surfaces of the tube. As shown in figure 2, if a 
beam of unit irradiance enters the tube wall from the 
inner surface, the total transmitted irra 
be represented by the following infinite 
appendix A): 
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(a) Region 1 — rays penetrate inner surface. 
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(b) Region 2 — rays totally reflected from inner surface. 



(c) Region 3 — rays never reach the inner surface. 

Figure 1. Light beams incident on a tube. For clarity, some 
transmitted beams are not shown. 



Figure 2. Typical beam path with multiple reflections and 
transmissions. 


reflectance at inner surface 

r Q reflectance at outer surface 

t 2 transmission at inner surface, 1 — ri 

t 0 transmission at outer surface, 1 — r a 

(These variables are functions of the polarization of 
incident radiation. Though not explicitly shown, the 
contributions of parallel and perpendicular polariza- 
tions have been included in the subsequent computa- 
tions.) Similarly, the total inside reflected irradiance 
Ri is given by the series 


R l = r i + t 2 x e 2 r 0 (l + e 2 r l r 0 + e A r 2 r 2 0 + ...) 


= n + 




(5) 


1 — e 2 r l r 0 

The total absorbed irradiance in the tube wall is thus 


1 -Ri-T (6) 

The absorbed irradiance of equation (6) applies 
to one infinite series of reflections between inner and 
outer tube walls as shown in figure 2. Actually, there 
are an infinite number of these infinite series. (Each 
beam from the original series that reenters the tube 
interior creates another such series with exactly the 


i 
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same geometry but with less irradiance.) To continue 
the analysis, let all the beams from the original 
series that reenter the tube interior be considered as 
one beam with irradiance /?, . This collective beam 
reflects around the tube inner surface, and at each 
reflection its irradiance is reduced by R j. Thus, 
the absorbed irradiance from the infinite number of 
infinite series is 


(1 -Rj - T) (l + R, + Rj + R'f + .-) 

1 -Rj-T 

1 ~Ri 


Recall that the analysis began with a beam of 
unit irradiance entering the inner surface of the tube. 
This beam must be related to an external irradiance 
if the absorption formulas above are to be related 
to external irradiance. A beam of unit irradiance 
incident on the exterior of a tube will transmit a 
fraction T to the interior. (The value T was derived 
above for a beam coming out of the tube but its value 
is independent of direction of travel.) Therefore, if 
the original beam incident on the inner surface has 
a value of T instead of unity, it may be regarded 
as the result of unit beam irradiance external to the 
tube. Accordingly, equation (7) must be multiplied 
by T. The absorbed irradiance for an infinite number 
of infinite series of reflections becomes 


T - 


T 2 

i - R, 


( 8 ) 


of incidence and 82 is the angle of refraction (see 
appendix B), total absorbed power per unit surface 
area becomes 


a(\, 9) 


1 -Ro 


T 2 \ cos 9 1 
1 — R{ ) cos 02 


( 11 ) 


Equation (11) applies to region 1 (fig. 1(a)), which 
results in rays penetrating the inner surface of the 
tube wall. In region 2 (fig. 1(b)), where there is 
internal reflection at the tube inner surface, the total 
absorbed irradiance is 

a(X,9) = (l-R 0 ) (12) 

cos 9 0 

In region 3 (fig. 1(c)), the total absorbed irradiance 
is given by 


a(A,0) 


(1 - e)(l — r„) cos 6 \ 
(1 - er 0 ) cos 9o 


(13) 


The calculation of e above requires a knowledge 
of absorption length in the glass d. Given R 2 , r, n, 
and 9\, as illustrated in figure 3, d can be shown to 
be 

d — i?2 cos 62 — \[ (Ro cos 0 2 ) 2 — 2 R 2 t + r 2 (14) 
where 

R 2 outside radius of tube 

r wall thickness 

0 2 sin -1 [sin(0]/ri)] 


With the addition to equation (8) of the irradi- 
ance absorbed from the series of reflections generated 
by unit irradiance incident from outside the wall, the 
absorbed irradiance a(A,$) becomes 


n index of refraction 

For region 3 of figure 1(c)), 

d = 2R 2 cos 9 2 (15) 


n(\.0) = (1 -R (i -T) + 



= 1 -Ro 


T 2 

1 - R, 


T 2 \ 

1 -Rij 


(9) 


where A refers to the spectral wavelength and the 
total outside reflectance R (> is 


Ro = r (l 


t 2 e 2 r 
*o c l l 


1 




( 10 ) 


Because incident irradiance decreases by the factor 
cos 9\j cos 80 inside the glass, where 8 \ is the angle 


Calculations 

The absorbed irradiance in the tube wall is a func- 
tion of the path length in the tube wall, the index of 
refraction of glass, the angle of incidence of the ir- 
radiance, the absorption coefficient of the glass, and 
the incident solar intensity. Intensity, index of re- 
fraction, and absorption coefficient are spectral vari- 
ables. Thus, absorbed irradiance must be integrated 
over the wavelengths at which absorption and solar 
irradiance are significant (i.e. , 155 to 4900 nm). Path 
length and reflection at tube boundaries are functions 
of the angle of incidence at the tube surface and the 
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Figure 3. Beam path length in glass wall. 

refractive index. Hence, absorbed irradiance must 
also be integrated over many angles of incidence (be- 
tween 0° and 90°) at each wavelength. Essentially, 
the following integral was evaluated: 


2 fX r? 

It = -J j S(\)a(X, e ) dO d\ W/m 2 (16) 


where 


5(A) 

solar spectral intensity distribution, 
W/m 2 -nm 

a( A, 6) 

spectral absorption (eq. (11), (12), 
or (13), depending on the value of 
Oi) 

Solar spectral intensity data were obtained from 
reference 6, and refractive index data were obtained 


from reference 7. Both references present data at 
wavelengths with nonuniform intervals. Transmis- 
sion data were obtained from reference 7 by digitiz- 
ing graphical data presented at irregular intervals. 
Intermediate values of these three parameters were 
obtained from a second-order interpolation of data 
at the next smaller and the next two larger wave- 
lengths. The absorbed irradiance was calculated for 
a range of tube wall thicknesses for two tube radii of 
50 and 5 cm. 

Results and Discussion 

Table I shows the refractive index of fused silica, 
the total reflectance (both surfaces) of fused silica, 


the absorption coefficient of fused silica, and the 
solar spectral intensity at various wavelengths. In 
the spectral region 155 to 220 nm, total reflectance is 
unusually large because the refractive index is large. 
The absorption coefficient has a very small value 
where the solar spectral intensity reaches significant 
values. Thus, absorbed irradiance in the ultraviolet 
spectrum is negligible. 

In the spectral region 220 to 2780 nm, values 
of the absorption coefficient are zero (within the 
accuracy of measurement). Practically no irradiance 
is absorbed in this spectral region, so no data are 
presented in table I. 

From 2780 to 4900 nm, by far the major portion 
of irradiance is absorbed. The particular glass chosen 
for these calculations minimizes absorbed irradiance 
because a large absorption band at 2700 nm has been 
eliminated by the removal of OH radicals from the 
glass. 

Beyond 4900 nm, 7 W/m 2 of irradiance are in the 
solar spectrum and the absorption coefficient is very 
large but unknown. It is estimated, therefore, that 
1 W/m 2 is reflected and 6 W/m 2 are absorbed. Thus, 
6 W/m 2 have been added to each calculated value of 
absorbed irradiance for the tabulated thickness. 

Values of absorbed irradiance for various tube di- 
mensions are listed in table II and plotted in figure 4. 
The values for large tube wall thicknesses can be used 
to approximate power absorption in solid rods. The 
absorbed irradiances expressed by the two curves at 
the same thickness show the effect of the radius of 
curvature to be almost negligible. 

Since equations (11), (12), and (13) for absorbed 
irradiance include the independent variables e, r 
and r 0 , it is surprising that the absorbed irradiance 
does not vary with rj! Physically, this means that 
the tube absorbs as if the interior surface were not 
there. Only rays that leave the outer surface affect 
absorbed irradiance. The inner surface serves only 
to distribute absorption around the tube. 

The calculations and these statements are made 
with the assumption that the tube interior contains 
no material that absorbs infrared radiation. The 
presence of infrared-radiation absorbing material in 
the tube interior would reduce absorbed irradiance 
in the walls of the tube. 

The calculated values of absorbed irradiance in an 
empty tube are minimum values for two significant 
reasons: (1) only rays in a plane perpendicular to the 
tube axis have been considered, and (2) absorption 
of the glass increases with temperature (ref. 8). In 
the most general case, incident rays will not be 
perpendicular to the tube axis, and the absorptive 
path length in the glass will be increased. This will 
increase the absorbed irradiance. 
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Table I. Spectral Properties of Fused Silica and Solar Spectral Intensity 


Wavelength, 

nm 

Refractive 

index 

Total reflectance 
(two surfaces) 

Absorption 

coefficient 

Solar spectral 
intensity 

155 

1.641361 

0,111353 

24.042201 

0.000120 

160 

1.628813 

.108240 

10.440855 

.000230 

165 

1.616886 

.105289 

5,607232 

.000403 

170 

1.605579 

.102498 

2.634072 

.000630 

175 

1.594894 

.099867 

1.280840 

.000835 

180 

1.584830 

.097397 

.600185 

.001250 

185 

1.575387 

.095085 

.236414 

.001164 

190 

1.566565 

.092931 

.046188 

.002710 

220 

1.528498 

.083718 

.002016 

.057500 

2780 

1.424320 

.059447 

.003927 

.039800 

2785 

1.424213 

.059423 

.008362 

.039600 

2790 

1.424106 

.059399 

.012766 

.039400 

2795 

1.423999 

.059375 

.017140 

.039200 

2800 

1.423891 

.059351 

.021484 

.039000 

2850 

1.422801 

.059107 

.064847 

.037000 

2900 

1.421684 

.058857 

.104285 

.035000 

2950 

1.420542 

.058601 

.133458 

.033125 

3000 

1.419370 

.058340 

.151974 

.031000 

3050 

1.418169 

.058072 

.158920 

.028300 

3100 

1.416940 

.057798 

.141714 

.026000 

3150 

1.415686 

.057518 

.127101 

.024300 

3200 

1.414400 

.057232 

.115086 

.022600 

3250 

1.413081 

.056939 

.070738 

.020800 

3300 

1.411730 

.056640 

.057293 

.019200 

3350 

1.410345 

.056333 

.076124 

.017825 

3400 

1.408930 

.056020 

.121736 

.016600 

3450 

1.407485 

.055701 

.207387 

.015488 

3500 

1.406010 

.055375 

.298713 

.014600 

3550 

1.404505 

.055044 

.499959 

.014063 

3600 

1.402970 

.054706 

.902126 

.013500 

3650 

1.401405 

.054362 

1.266029 

.012900 

3700 

1.399810 

.054012 

1.589068 

.012300 

3750 

1.398185 

.053657 

1.809609 

.011650 

3800 

1.396530 

.053295 

1.878325 

.011100 

3850 

1.394845 

.052927 

1.900677 

.010700 

3900 

1.393130 

.052554 

1.915146 

.010300 

3950 

1.391385 

.052175 

2.057062 

.009900 

4000 

1.389610 

.051790 

2.336959 

.009500 

4050 

1.410607 

.056391 

2.752873 

.009113 

4100 

1.428891 

.060475 

3.277175 

.008700 

4150 

1.444462 

.064004 

4.005107 

.008225 

4200 

1.457321 

.066952 

4.798062 

.007800 

4250 

1.467468 

.069297 

5.664125 

.007438 

4300 

1.474901 

.071026 

6.770718 

.007100 

4350 

1.479623 

.072128 

7.797080 

.006798 

4400 

1.481632 

.072598 

8.756014 

.006500 

4450 

1.480928 

.072434 

8.966105 

.006209 

4500 

1.477512 

.071635 

8.532039 

.005920 

4550 

1.471383 

.070207 

7.979737 

.005625 

4600 

1.462541 

.068156 

7.646193 

.005350 

4650 

1.450988 

.065496 

8.391921 

.005093 

4700 

1.436721 

.062244 

9.640547 

.004860 

4750 

1.419742 

.058423 

11.822024 

.004661 

4800 

1.400051 

.054065 

17.219726 

.004470 

4850 

1.377647 

.049214 

21.249766 

.004285 

4900 

1.352530 

.043925 

33.691100 

.004110 









Figure 4. Variation of absorbed irradiance per unit surface area of fused silica (OH free) glasi 


The temperature reached by the tube wall will de- 
pend on the spectral emission coefficient. That is, in 
the absence of any other cooling, the tube will reach 
a temperature at which radiated power will equal ab- 
sorbed irradiance and, since the absorbed irradiance 
will be multiplied by the solar concentration factor C 
at the tube surface, 

CI t = CSa^ ecrT? (17) 

where 

e effective emissivity 

S effective solar power 

a Stefan-Boltzmann constant, 

5.7 x 1(T 8 W/m 2 — K 4 

a effective absorptance 

Tt tube temperature 

Thus, 


The evaluation of equation (18) is complicated 
because e is also a function of temperature. In con- 
trast to the surface emission of opaque materials, the 
tube radiates energy throughout its volume. Gardon 
(ref. 9) has shown that the emissivity of volumetric 
radiators depends on the product of thickness and ab- 
sorption coefficient and has evaluated the total hemi- 
spherical emissivity of ordinary plate glass for several 
thicknesses. If the emissivity of fused silica approx- 
imates the emissivity of plate glass, the operating 
temperature of the tube can be estimated. 

For purely radial rays, solar concentration C can- 
not exceed 215 (see appendix C). For this concen- 
tration, an iterative solution to equation (18) with 
the data of Gardon for a 0.2-cm-thick tube shows the 
temperature to be 519 K. Conversely, if the tube does 
not exceed the maximum continuous operating tem- 
perature of 1225 K, the solar concentration cannot 
exceed 4450. (In this case rays must have longitudi- 
nal as well as radial components.) 

The above temperature of 1225 K is conservative 
and the concentration is liberal because the increase 
of a above 700 K (ref. 8) has not been included. This 


is mainly a surface effect that is absent in the pres- 
ence of oxygen. In a vacuum, the effect will cause 
surface heating that will raise the tube temperature 
by about 25 percent and lower the permissible so- 
lar concentration. Conductive cooling may occur to 
some degree if the lasing medium contacts the tube 
walls, or it may be introduced intentionally by other 
means. Regardless, conductive cooling will reduce 
heating to some degree but will create thermal gra- 
dients in the walls of the tube that, if large enough, 
can damage the tube. Wall temperatures can be re- 
duced by reducing the infrared (IR) irradiation of the 
tube in several ways: (1) the lasing medium could ab- 
sorb part of the IR radiation; (2) the tube could be 
coated with an IR reflector; or (3) the IR radiation 
could be absorbed or not reflected before it reaches 
the tube. 

Conclusions 

The formulas in this report provide a means for 
calculating absorbed power (irradiance) in any glass 
tube and estimating the operating temperature of the 
tube. The absorbed irradiance varies primarily with 
tube thickness and is independent of the reflectivity 
of the inside surface. Variation of absorbed irradi- 
ance with tube diameter is insignificant. 

In the application of the formulas to tubes made 
of hydroxyl-free fused silica it was found that, with 
radiative cooling only, the tube would operate at ap- 
proximately 519 K with purely radial irradiation and 
could reach its maximum operating temperature at 
a solar concentration of approximately 4450. Irra- 
diance is absorbed almost entirely at infrared (IR) 
wavelengths. Since this form of fused silica absorbs a 
minimum of IR radiation, other forms of fused silica 
can be expected to get hotter or tolerate less solar 
concentration under the same conditions. 

Tube heating can be lessened by (1) partial ab- 
sorption of IR radiation by the lasing medium, (2) 
coating the tube with an IR reflector, (3) removing 
IR radiation from the sunlight before incidence, and 
(4) using a different material that does not absorb 
significantly. Conductive cooling can lessen the heat- 
ing but will create thermal gradients in the glass that, 
if large enough, can damage the tube. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
February 25, 1987 
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Appendix A 

Transmitted Intensity 

Assume a light ray of unit intensity enters the inner surface of a tube from the interior. At the inner surface 
a portion t j of the ray intensity is transmitted into the tube wall. That portion is attenuated by the factor 
e before reaching the outer surface. At the outer surface part of the reduced intensity t t e is reflected back 
toward the inner surface and part is transmitted outside the tube. The part transmitted outside the tube is 
equal to tfi 0 e and is the first part of the transmitted intensity. The portion reflected travels back to the inner 
surface beginning with intensity t 2 er 0 and arrives at the inner surface with intensity tje 2 r c . After reflection 
at the inner surface the intensity becomes <je 2 r 0 rj and proceeds, again, to the outer surface. The intensity 
arriving at the outer surface is reduced again by e to fie^r 0 ri, and tfi 0 e 3 r 0 r^ is transmitted outside the tube. 
The transmitted intensity becomes 


T = tfi 0 e + tit 0 e 3 r 0 ri +•■■ 

= tit 0 e ^1 + e 2 r 0 r 2 + •• (Al) 

At any point on the outer surface where a transmitted beam emerges there is also a reflection which produces 
the next transmitted beam. Each succeeding transmitted beam is reduced in intensity by e 2 r 0 r^. Hence, the 
series above for the total transmitted intensity can be continued ad infinitum. Since 

1 + e 2 r 0 ri + e 4 r 2 r? + e 6 r^rf + • • • = (e 2 r 0 r; < 1) 

l — e*r 0 r ? - 


then 


T = 


t{t 0 e 

1 - e 2 r 0 ri 


(A2) 
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Appendix B 

Variation of Radiance With Angle of 
Incidence 

Let #1 be the angle of incidence on a solid surface 
{6 1 is measured from a perpendicular to the surface). 
At normal incidence (6\ = 0 ), any infinitesimal spot 
of light on the surface covers a minimum area gL 4 . As 
6\ increases, the area elongates in the direction of the 
plane of incidence. This elongation causes the area to 


increase by the factor 1 / cos 9 \ , and spot irradiance 
decreases by cos 0\. 

Inside the solid material, the angle of refraction 
62 is smaller than 6\ (Snell’s law) and decreases 
the one-dimensional elongation caused by 6\. The 
area perpendicular to the direction of beam propa- 
gation in the material decreases by the factor cos 62, 
so that the beam area inside the solid becomes 
dA cos 62 / cos 6\ . Since irradiance is inversely pro- 
portional to the area of the beam, irradiance de- 
creases by cos 6 1 / cos 62- 
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Appendix C 

Maximum Solar Concentration in One Dimension 

If sunlight is incident on the frustum of a reflective right circular cone through the annular area 7r — rfj 

in its base and the apex angle of the cone is 90°, then the rays will be reflected to the smallest possible region, 
which is the frustum of a much smaller cone centered on the axis of the reflective cone. (See fig. Cl.) The 
radius of the small cone at its base is r 2 <a and at the other end is rja (a is the half-angle of the solar subtense 
angle, 4.6543 x 10~ 3 rad). The length of the small cone is T2 — ri from the geometry. The area of the frustum 
of the small cone is given by 


7r(ria + r 20 t)\J [r% — ri) 2 + {r^cx. — rice) 2 (Cl) 

Maximum solar concentration C is thus expressed as the ratio of the annular area to the frustum area of the 
small cone: 


C = 


- ( J2 

„2\ 




n { r 2 



1 

1 

Cv 


7r(ric* + r 20 i)yj (j 2 — ri) 2 + (r 2 <* — ria) 2 a\/l + a 2 a 


ss 215 


(C2) 



Figure Cl. Concentration of sunlight by a conical reflector into a conical focal region. 
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